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Abstract Pulsar glitches, i.e. the sudden spin-ups of pulsars, have been detected for 
most pulsars that we known. The mechanism giving rise to this kind of phenomenon is 
uncertain, although a large data set has been built. In the framework of star-quake model, 
based on lBavm & Pinesl dl97ll ). the glitch-sizes (the relative increases of spin-frequencies 
during glitches) Af2/f2 depend on the released energies during glitches, with less released 
energies corresponding to smaller glitch sizes. On the other hand, as one of dark matter 
candidates, our Galaxy might be filled with the so called strange nuggets (SNs) which 
are the relics from the early Universe. In this case the collisions between pulsars and 
SNs is inevitable, and these collisions would lead to glitches when enough elastic energy 
has been accumulated during the spin-down process. The SNs-triggered glitches could 
release less energy, because the accumulated elastic energy would be less than that in 
the scenario of glitches without SNs. Therefore, if a pulsar is hit frequently by SNs, it 
would tend to have more small size glitches, whose values of ACl/Cl are smaller than 
that in standard star-quake model (with larger amounts of released energy). Based on the 
assumption that in our Galaxy the distribution of SNs is similar to that of dark matter, as 
well as on the glitch data in ATNF Pulsar Catalogue and lodrell Bank glitch table, we find 
that in our Galaxy the incidences of small size glitches exhibit tendencies consistent with 
the collision rates of pulsars and strange nuggets. Further test of this scenario is expected 
by detecting more small glitches (e.g., by the Square Kilometre Array). 


1 INTRODUCTION 

Glitches are one type of pulsar timing irregularity, observed as discrete changes on the pulsar rotation 
rate. Although a large data set for glitches has been built, the physical picture behind them still remains 
to be well understood. Generally, there are two main models for the origin of glitches. The first one is the 
star-quake model which regards a_glitch a s a s tar-quake of a spinning down pulsar resulting from the re¬ 
arrangements of its crust (iBavm & PinesLl 197 it) . The second one regards glitches_as the result of a rapid 
transfer of angular momentum between the inner superfluid and the outer crust d Anderson & Itohll 19751) . 
Although the star-quake mod el is now less popular which falls to giv e a reasonable explanation for the 
glitch activity of Vela pulsar (Bav m & Pineslfl97lUAlpar et al.Ul996h . the rearranges o f the crust could 
be th e trigger for the transfer of angular momentum described in the second model (Esp inoza et all 
1201 ll) . What is certain to us is that, glitch behaviors reflect the interior structure of pulsars, and the data 
of glitches can give important information on the dense matter inside pulsars. Some correlations between 
the glitch activities and the characteristic age of pulsars have been found, e.g., younger pulsars would 
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have larger glitch spin-up rates O and glitch-rates N g (iEspinoza et al.Li2071 ). The data of glitches make 
us to believe that glitches are caused by the changes of interior structures. Nevertheless, can glitches 
related to some exterior trigger? 

Another topic, seemly having no relation with glitch phenomenon, is about the relics of the early 
Universe. In the Standard Model of cosmology, during the expanding history, our Universe underwent 
several kinds of phase transitions. The QCD (quantum chromodynamics) phase transition, taking place 
when the age of the Universe was about 1 millisecond and the temperature was about 200 MeV, would 
leave out a huge amo unt of relics, if the phase transition was of the first order dWittenl 1 1984 iMadseni 
Il999h . iwittenl (il984l) suggested that most of the deconfined quarks would be concentrated in the dense 
and invisible strange nuggets (SNs). In fact, the formation and evolution of SNs depend on the physical 
properties of hot quark-gluon plasma and the state of dense matter, so the exact baryon number A 
of each SN (or the mass-spectrum of SNs) is hard for us to det e rmine . Although some calculation s 
were made, e.g. iBhattachariee et all dl993l> : iBhattacharvva et all (l2000l) : iLugones & Horvathl d20Q4l) . 
the results were model dependent. On the other hand, the properties of SNs could be constrained by 
astrophysics, if they do exist in the Universe. 

As proposed by Witten, SNs could be a kind of dark matter candidate dWittenl 119841) . Although 
being baryonic matter, SNs could behave like dark matter because of the very low charge-to-mass ra¬ 
tio, which comes from the fact that SNs are made of nearly equal numbers of it, d and s quarks. As 
will be shown in this paper, the mass of each SN could be as high as ~ 10' g (with baryon num¬ 
ber A ~ 10 :! 1 and size ~ 10 -2 cm). Whether such macroscopic and baryonic particles are com¬ 
patible with observations, such as the formations of galaxies and stars, should be studied in some 
detail, but unfortunately even the numerical simulations in the framework of conventional micro¬ 
scopic dark_matter particles have not reached consensus about the cosmic structure formation (e.g. 
see lGao & Tlieuns ( 200 7)). A similar macroscopic dark matter candidate are the so called massive com¬ 
pact halo objects (MACHOs) dPaczvnskil Il986l) . and some constraints from observations have been 
found. Very heavy MACHOs with masses above 10~ 7 M^ are excluded to be primary dark matter par¬ 
ticles dTisserand et all 120071: iMetcalf & SilkL 1 20071) . but dark matter could be made out of baryons 
and antibaryons in the form of com pact stellar type objects or primordial black holes with suffi- 
ci ently wide mass spectrum (Dolgov et all l2009h . The latest constrains on MACHOs could be found 
in Monrov-Rodriguez & A Hern 2014 ). which put the upper bound for the mass to be 5 Mq. The stud¬ 
ies about MACHOs could infer that, whether macroscopic dark matter particles with masses below 
10~ 7 M© have any conflict with either observations or simulations remains to be seen, and the related 
problem about how the SNs affect the structure formation is worth exploring in the future. 

SNs as dark matter should not lead to conflicts with observations, and on the other hand the observa¬ 
tions could give constraints on the properties of SNs. lMadsen & Riisagerid 19851) found that if dark matter 
is completely composed of SNs, the consistency between the prediction of Big Bang Nucleosynthesis 
(BBN) and the observational result for the helium abundance would not be affected if the radius of 
SNs exceed ~ 10 -5 cm, corresponding to the baryon number A of each SN with A > 10 24 (assum¬ 
ing that the matter inside SN has the density comparable to that of nuclear matter density). SNs in the 
early Universe would affect the structure formation. On the other hand, SNs could help us to understand 
some astrophysical phenomenons. The formations of high-redshift (z ~ 6) supermassive black holes 
(10 9 Mq) dFan et all 120031) could be the results of SNs dLai & Xull2010h . and in this picture the baryon 
number of each SN should be smaller than about 10 35 . Another astrophysical consequence of SNs in our 
present Universe is that, i f they pass our Earth, some seismic signals would be detected (And erson et all 
[20031 : [Herrin et al.ll2006l) . More detailed studies about the effects of strange nuggets on some important 
processes, such as the formation of early galaxies and black holes, are expected. 

In this paper, we propose another astrophysical consequence, which is related to the pulsar glitches. 
If SNs formed in the early Universe and survived in the cosmic evolution, they would fill our Milky 
Way, and they would inevitably be accreted by stars. With strong gravitational field, the accretion of 
SNs by pulsars would not be rar e events. In fact, the accretion of dark m atter particles by sun as well 
as by pulsars have been studied dPress & Snergel. 1985: Kouvaris , 2008 ), and som e interesting conse¬ 
quences have been discussed, i.e ( Goldman & Nussinov , 1989t Huang et al., 12014 Assuming that in 
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our Galaxy the distribution of SNs is similar to that of dark matter, we can then know the accretion rates 
of SNs by pulsars. Whether such accretion would affect the properties of glitches depends both on the 
states of matter inside pulsars and the glitch mechani sm. In thi s pape r, we assume that pulsars are solid 
strange stars and apply the star-quake model dBavm & PinesLIl97ll) as the glitch mechanism. In fact, 
the star-quake glitch m odels of solid strange stars have been studied (Jzhou et all 120041 : iPeng & XuL 
2008; IZhou et all 120141) . which show that glitch behaviors and the energy released are consistent with 
observations. Moreover, in star-quake model, the recovery process after a glitch could be explained as 
the damped oscillations (IZhou et al1120041) ._ 

In the framework of star-quake model dBavm & Pinesl 1T97 ll) . the relation between glitch sizes (the 
relative increases of spin-frequencies during glitches) Afi/fl and the released energy during glitches 
can be derived straightforward, which shows that less amount of released energy during glitches would 
correspond to smaller glitch sizes. Pulsar structure would be affected by the collisions between the pulsar 
and SNs, which could help the pulsar to release its elastic energy that accumulated during the spin-down 
process, probably manifesting as glitches detected. In this scenario, the amount of energy released in 
this glitch E re would be smaller than that in the scenario without SNs, as the energy released by the SN 
itself is negligible. Therefore, for a particular pulsar, the incidence of small size glitches should be the 
same as its accretion rate of SNs. Here the small size glitches are relatively small ones, whose values of 
Afi/fl are smaller than that in the standard star-quake scenario without SNs. 

In this paper, we demonstrate the relation between accretion of SNs and the incidences of small size 
glitches in the star-quake model. Based on the assumption that in our Galaxy the distribution of SNs is 
similar to that of dark matter, as well as on the glitch data in ATNF Pulsar Catalogue and Jodrell Bank 
glitch table, we find that in our Galaxy the incidences of small glitches exhibit a tendency consistent 
with the collision rates of pulsars and strange nuggets. It should be noted that, although in this paper we 
assume pulsars are totally in solid state, our conclusions would not change quantitatively if pulsars only 
have solid crusts. 

This paper is arranged as follows. In §2 we will derive the relation between the glitch sizes and the 
released energies during a glitch, in the star-quake of solid strange star model. Based on this, the relation 
between accretion of SNs and the small sized glitches will be demonstrated in §3. In §4 we will get the 
expected accretion rates of SNs by pulsars and compare that with the observed event rates of small size 
glitches. Conclusions and discussions are made in §5. 


2 GLITCHES IN STAR-QUAKE MODEL 

Following iBavm & Pinesl d 197 lb . we will derive the released energies in terms of the glitch sizes (the 
relative increases of spin-frequencies during glitches) AH/H. In addition, we do the calculation in the 
framework of solid strange star model, which means that all parts of the stars are in solid state, but the 
results would not change quantitatively in the model of neutron stars with solid crusts. When a solid star 
spins down, strain energy develops until the stresses reach a critical value. At that moment, a star-quake 
occurs and the stresses are relieved and its rotation rate is suddenly increased by conservation of angular 
momentum (similar calculations can be found in IZhou et al ( (20041) an d IZhou et al.lJ2014k 
Total energy of a rotating star with mass M, radius R and angular momentum L is 


E — E 0 + — + Ae 2 + B(e — eg) 2 , 


where E 0 is the total energy in the non-rotating case, and 


I-Io 

h 


( 1 ) 

( 2 ) 


is the oblateness of the star with moment of inertia /, which reduces to / (J = 0.4 A/ li 2 in the non¬ 
rotating case, and eo is the initial oblateness. The coefficient measure the departure of gravitational 
energy relative to the non-rotating case is A = ^GM 2 /R, and the coefficient measure the strain 
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energy is B = i^V C/u where C44 is the shear modulus and V = IWi’ 3 /.') is the whole volume of the 
star dBavm & PinesLI 197 ll) . The angular momentum is conserved during a glitch, so the sudden chance 
in oblateness Ae is directly related to the relative jump in angular momentum fl 


Ae = e-i - e+» = - 


A D_|_j — ri¬ 


ft 


+i 


fl 


+i 


(3) 


where the spin-up rate mentioned above Aft/ft is just Aftj/ft + ,- during i-th glitch, and the subscripts 
“+*” and i” denote quantities right before and right after i-th glitch. The initial oblateness before v-th 
glitch e oi is dBavm & PinesL H97lh 


/ 0 ftL b 

e+i - 4 (A + B) + ATb £ oi ' (4) 

The evolution of e as time is shown in lZhou et al.l ( 2004 1. 

For simplicity we assume that when a star-quake occurs, the entire stresses are relieved and the 
released energy is E re . The energy conservation gives 

2 ^+tfl+i + Ae 2 +i + B(e + i — eoi ) 2 = -/-jftlj + Ae 2 i + E re , (5) 

then it is easy to get 

Af li _ E re )j, _ 2 E re (A + B ) 1 

fV~~2 ~b^ c ~ bT 0 n 2 t - n 2 +i ’ 

where a c is the critical mean stress of the star, which corresponds to the occurrence of a glitch, and 
[i = 2B/V is the mean shear modulus of the star. 

The value of parameter B for solid strange stars is uncertain now, because the structure and state of 
quark matter inside pulsars are still unsolved problems. The shear modulus C44 for a bee lattice of nuclei 
of number density n, charge Ze a nd lattice constant a, interacting via an unscreened Coulomb interac¬ 
tion, could be written as dBavm & PinesLi~j~97~ 1) C44 ~ 0AZ 2 e 2 n/a ~ IO 30 erg/cm 3 . Based on the fact 
that the strong interaction is about 2 to 3 orders of magnitude stronger than electromagnetic interaction, 
we infer that the value of C44 of solid strange stars is about in the range 10 32 — 10 33 erg/cm 3 dXul 
l2003ll . then B = ^VC 44 ~ 10 52 erg ■ (C 44 /IO 33 erg/cm 3 ) (R/10 km) 3 . Comparing with A ~ 
6 x 10 52 erg (M /1.4 M Q ) 2 ■ (10 km /R), we can approximate that (A + B)/B ~ 10, then 


Aflj ~ E ie (A + B) 1 

f 2+i BI 0 n 0i |f 2 |Af 


~ 10 


-6 


E r , 


10 34 erg 


10 45 


g • cm 


10 




(7) 


where P = 2tt/ f > is the spin period, and Af is the time interval between two successive glitches. The 
result is consistent with Fig.2 in IZhou et al.l d2014l) which shows the upper limit of the released energy 
during glitches (for bulk invariable case). 


3 STRANGE NUGGETS AND GLITCHES 

3.1 Cosmological QCD phase transition and strange nuggets (SNs) 

Cosmological QCD (quantum chromodynamics) phase transition, taking place when the age of the 
Universe was about 10 6 seconds and the temperature was about 200 MeV, could have interesting 
astrophysical consequences. If it was of the first order, a huge amount of strange nuggets (SNs) would 
form and survive the cooling of the Universe, and SNs were proposed as dark matter candidates ( WittenL 
119841) although they are baryonic, due to their extremely low charge-to-mass ratio (a detailed review 
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of SNs could be found in iMadsenl 1 1999 1). SNs are composed of up, down and strange quarks, with 
densities higher than nuclear matter density. The state of SNs is uncertain, because this problem is 
related to the low-energy QCD theories, the same problem about the interior structure of pulsar-like 
compact stars. 

SNs with small baryon number A would have been evaporated in the early Universe, but the lower 
limit for the baryon number A of the stable SNs is difficult for us to derive. On the other hand, some 
constraints on size of SNs have been made combined with their astrophysical consequences. If the 
existence of SNs would not break the consistency between the prediction of Big Bang Nucleosynthesis 
and observed result for the helium abundance, the baryon number of each SN should be larger than 
about 10 24 (assuming that SNs have a uniform size) ( Madsen & Riisageil Il985l). Using some m ore 
accurate data, the lower limit of baryon number has been improved to be ~ 10 25 (Lai & Xu blllOl) . If 
the formations of high-redshift (z ~ 6) supermassive black holes (10 9 Mq) dFan et all 12003 ) were the 
results of SNs, then the baryon number of each SN should be smaller than about 10 35 ( Lai & Xu[l20lO ). 

Then what could SNs do in our present Galaxy? As shown in IWitteri i 19841) . the mass fraction 
of SNs could be as large as that of dark matter. If SNs are stable and survive in present Universe, they 
would fill our Galaxy, with the total mass comparable with that of dark matter. Circulating in the Galaxy, 
they would be accreted by stellar objects, especially pulsars whose gravitational fields are strong, and 
the accretion rate by a particular pulsar is then related to its position in the Galaxy. 

Then we will derive the accretion rates of SNs by pulsars. Assuming that SNs have the same distri¬ 
bution as that of dark matter, and the fraction of total mass of SNs and compared to the total mass of dark 
matter is r /, whose value is chosen to be 0.1 in the following calculations. For the dark m_atter density 
distribution and the function of distance from Galaxy center r, we choose the NFW form ((Navarro et all 
119971) 

Ps 


Pdm(r) = 


(*) 1+ fe)' 


( 8 ) 


where p s = 0.26GeV cm' 3 and r s = 20 kpc. The number density of SNs is n x = rj ■ p& m /m x , 
where m x is the mass of each SN. For simplicity, we assume that all of SNs have the same mass, with 
m x = A ■ mb, where A is the baryon number of each SN and mi, ~ 1 Ge V is the mass of baryon. The 
circular velocity v of each SN is then 


j(r) = 


GM(r) 


(9) 


where M (r) = f ( ' p{r' )47 t?’' 2 dr' is the total mass inside the sphere with radius r. To derive the accretion 
rates of SNs by pulsars, we can borrow the results for the accretion rates of dark matter particles by 
sta rs/pulsars. The accretion rate of SNs wi th number density n x by a pulsar with mass A1 and radius R 
is dPress & SpergelL[T985UKouvarisil2008l) 


T = 47 xn x 


3 \ 3/2 2 GMR 1 

2t tv 2 J 1 - 2 GM/R 3 * 


GO) 


where the factor (1 — 2GM/R) 1 is the modification factor taking into account the general relativity 
effect. 


3.2 SNs and glitches 

The sample of all of the observed glitches also shows that, the glitch sizes Afl/fl spread between 
1CU 11 and 10j 4 with a bimodal distribution which peaks at approximately 10~ 9 and 1CU 6 , respec¬ 
tively (lEspinoza et all 1201 1 1:1 Yu et al.L 12013b . The glitch sizes have not been well explained, but Eq.® 
tells us that the glitch sizes are related to the released energies during glitches. When the stresses reach a 
critical value, the star releases strain energies E le , which are stored in the star during the spinning down 
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process due to the rigidity of the star. The released energy E le in a glitch can be inferred from observa¬ 
tions, if it is large enough for us to detect. For example, a constraint on the X-ray flux enhancement of 
Vela 35 days after a glitch with Afl/fl ~ 10 -6 have been made ( Helfand et al.Ll200H) . which infer that 
the upper limit of the energy released is about 10 36 erg. 

If a pulsar accretes an SN, its structure would be affected, as the SN would interact with particles 
inside the pulsar. The detail of such interaction is hard to derive, because it is related to both the interior 
states of pulsars and SNs. For an estimation we take the parameters used bv IZhou et ali ( 20041 . where 
glitches occur as the consequence of the star-quake of solid strange stars. When the stress reaches a 
critical value tx c , (e — cq) ~ 10 -3 , then the strain energy density p s train c °uld be approximated as 


Pstrain — C44 (fz 60) 

~ 10 -7 MeV • fm -3 


C44 


10 32 erg • cm -3 



( 11 ) 


If the lattice (quark or quark-cluster) density inside pulsar is ~ 0.3 fm~ 3 , then to break the lattice 
structure each lattice should obtain about 10 0 MeV. Just before hitting the surface of a pulsar, the 
energy of an SN is about (A/10 33 ) x 10 35 MeV, so if each lattice it goes through gains energy V, then 
the total number of lattices it would interact with is N ~ 10 41 • (A/10 33 ) • (10 -6 MeV /V). The cross 
section of the interaction between an SN and the particles inside the pulsar could be approximated as 
the geometric cross section of the SN, which is ~ 10 -4 cm 2 • (A/10 33 ) 2 / 3 , so the length of trajectory 
is ~ 10 km ■ (A/10 33 ) 1 / 3 • (10 ~ 6 MeV/V). 

Although the exact demonstration about the breaking criterion is unknown, the above estimation 
shows that it is possible for an SN with A > 10 33 to break the whole solid body (or crust) of the 
star, which would give rise to a glitch. The amount of energy releasing of such glitches triggered by 
SNs would be smaller than that of the normal glitches, as the latter occur when the elastic energy 
accumulated reaches a critical value. According to Eq©, in the star-quake model the glitch size Afi/fi 
is proportional to the amount of energy releasing E ie , so the sizes of SN-triggered glitches should be 
smaller than that of the normal glitches. 

The amounts of energy releasing E re in normal glitches are also difficult for us to quantify. It 
should be different from pulsar to pulsar, and could even be different between different glitch-events of 
one pulsar. The released energy is mainly from the accumulated elastic energy, and it could be inferred 
that the more massive pulsars should release more energy during glitches. This dependence of E le 
on the pulsar mass M could eliminate some uncertainties of both E ie and the moment of inertia Iq. 
Consequently, from Eqf7|) we can infer that, we could use typical values of both E le and Io to reduce 
the difference of AO/f2 between different pulsars. Because the uppe r limit of the energy released in on e 
glitch event of Vela (PSR B0833-45) pulsar is about 10 36 erg (Hel fand et all 1200ltlZhou et alll20l4 . 
it seems reasonable to set E le = 10 35 erg for normal glitches. 

Based on the above arguments, we could estimate the occur rates of SN-triggered glitches, which 
may be identified as those with AQ/fl < AQ/Cl(E le = 10 35 erg). On the other hand, we take Ea. dTOl) 
has the expected occur rates for SN-triggered glitches. Combined with the observational data, we could 
then compare these two rates for different pulsars. 

It is worth mentioning that, although the accretions of SNs by pulsars would lead to bursts in X- 
rays, it is not easy for them to be detectable. The gravitational energy released by accretion of an SN 
with A = 10 34 on a typical pulsar is about 10 3 ° erg. With degenerate electrons with Fermi energy 
~ 10 MeV, the heat conductivity of strange matter must be very high, which means that strange stars 
will be thermalized by any heat flow. The process of thermalization is complicated, but the time-scale 
T en h for releasing enhanced X-ray bursts resulting from the accretions could be estimated as r en h 
r ~ IQ 30 erg/ £hoi ~ 0.1 s • (10 31 erg s _1 /Lboi). where Lboi is the bolometric X-ray luminosity 
(see lYu & Xu d201 if) and references therein). The luminosity of an enhanced X-ray burst resulting from 
one accretion is then L e nh = IO 30 erg/r en h <C IO 30 erg/r ~ 10 31 erg s _1 , i.e., L en h is much smaller 
than the bolometric X-ray luminosity. Therefore, the accretions of SNs by some near and well-monitored 
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compact stars, such as RX J0720.4-3125 and RX J1856-3754, would lead to X-ray variability which are 
difficult to be detected, so the accretions would not bring conflict with observations. 


4 COMPARISON OF ACCRETION RATES OF SNS BY PULSARS AND GLITCH DATA 

As demonstrated in § [3] the collision between an SN and the pulsar could result to a small 
size glitches compared to that in the scenario without SNs. To compare with the glitch data, 
we choose 25 pulsars whose recorded numbers of glitches are larger than 5 ( Espinoz a et al 


[2011 ). The data for glitch sizes Afl/fl are from Jodrell Bank glitch-tables (jEspinoza et al. . 201 ll 


(http://www.jb.man.ac.uk/pulsar/glitches/html), and the rotation periods and the time derivation of pe¬ 
riods, as well as the locations of pulsars compared to the barycentre of the solar system (which 
are used to derive the distances of pulsars to the Milky Way center) are from ATNF Pulsar 
Catalogue dManchester et all 1200 5) (http://www.atnf.csiro.au/people/pulsar/psrcat/). 

The accretion rates of SNs by a pulsar (which are taken to be the expected event-rates of small 
glitches) with mass M = 1.4Mq and radius R = 10 km and located at a distance r from the Galaxy 
center are shown in Fig. 1, where solid, dashed and dash-dotted lines correspond to three values of baryon 
number A of each SN respectively. These lines denote the expected occur rates for small glitches with 
Afl/fl < Afl/f l(E le = 10 35 erg). The data points (blue stars) show the observed event-rates of small 
size glitches with Afl/fl < Af l/VL(E Te = 10 35 erg) and the distances from the Milky Way center 
of the corresponding pulsars. Here the event-rates are derived as the ratio of the number of glitches 
satisfying Afl/fl < Afl/fl(f? re = 10 35 erg) to the whole time interval from the first and the last 
observed glitches. The axes in Fig. 1 are scaled in logarithm, so the points corresponding to zero observed 
value are not shown in the figure, expect for J0537-6910. The observed value of J0537-6910 is shown 
on the horizontal axis although the value is zero, because it is the most distant one from the MW center 
among all these pulsars, and is consistent with our expectation. In the calculations, we assume that SNs 
constitute 10% of dark matter (?/ = 0.1). Although the mass ratio of SNs to dark matter // is unknown, it 
is proportional to A under a certain accretion rate. For instance, if // = 0.01, then the values of A should 
be one tenth of that in Fig.l. 

We can see that, except for the pulsars with zero event rate, the event-rates for glitches with 
Afl/fl < Afl/fl(E re = 10 35 erg) decrease with distances of pulsars from the Milky Way center, 
which show the same tendency as the accretion rates of SNs by pulsars. To demonstrate the correlation 
of glitch rates on the distances from the Galactic center, we use the Pearson correlation coefficient p, 
which is a measure of the linear dependence between two variables!]]. The result for data in Fig.l (ex¬ 
clusive of J0537-6910) is p ~ —0.5. Although the interpretation of a certain value of p will be different 
for different context and purposes without a criterion for all problems, the value —0.5 may imply the 
trend that, the more distant of a pulsar to the Galactic center, the less rates of small glitches it will have. 

The dispersion of data points might mainly result from two factors: (1) the fact that not all the SNs 
have the same baryon number A; and (2) the mass and radius of each pulsar would be different from the 
typical values. All of the 25 pulsars are listed in Table 1, including their spin-down ages, numbers of 
glitches (NGlt) and distances from the Milky Way center. The upper part of the table includes 9 pulsars 
who have glitch events with Afl/fl < Afl/f l(E le = 10 35 erg), and the lower part of the table includes 
16 pulsars who have no glitch event with Afl/fl < Att/tt(E re = 10 35 erg). 

Now we make some discussions about why some pulsars are out of our expectation, i.e. have no 
glitches satisfying Af1/fl < Afl/f l(E le = 10 35 erg). Comparing the data in Table[l] we can see that 
the spin-down ages of pulsars in the lower part of Table [j] are generally smaller than that in the upper 
part. As demonstrated before, we identify the SN-triggered glitches as the ones whose sizes satisfy 
Afl/fl < Afl/fl(i5 re = 10 35 erg), which means that we assume that all of the pulsars have the same 
amount of released energy during normal glitches. This assumption is certainly over simplified, and 
different pulsars should have different values of E re . The values of E ie of young pulsars would be 


1 For two variables x and y , p is defined as p = — = ~ =—( ' IJ ‘ = In general \p\ ^ 1, where p = 1 represents 

vEi(*r*) VEiGi-v) 

total positive correlation, p = 0 represents no correlation, and p = — 1 represents total negative correlation. 
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Fig. 1 The accretion rate of SNs, with baryon number 
A = 10 34 , A = 5 x 10 33 and A = 2 x 10 33 , by a pulsar 
located in the distance of r from the center of the Milky 
Way, are plotted in dashed, solid and dash-dotted lines 
respectively. In this paper, these lines denote the expect 
occur rates for pulsars with Afl/fl < Afl/fl(E re = 
10 35 erg). Data points: Observed rate per year of glitch- 
events with Afl/Q < Af l/Q(E ie = 10 35 erg), for 
pulsars whose recorded glitch-number are larger than 
5. The data point corresponding to J0537-6910 in the 
horizontal axis corresponds to zero observed value. 


larger than that of old pulsars, because the former ones are more active. The pulsars in the lower part of 
TableQ]might mean that we should impose some larger values of E le for the normal glitches. However, at 
present stage we are lack of a reliable description about the difference of E le between different pulsars, 
so we only use the criterion Afi/fi < A Cl/Cl(E Te = 10 35 erg) to identify the SN-triggered glitches. The 
result shown in Figure 1 is only the first attempt to reveal the correlation between the glitch behaviors 
and the location of pulsars inside the Milky Way. Whether the collision between SNs and pulsars could 
be a glitch trigger needs to be tested by further observations. 

5 CONCLUSIONS AND DISCUSSIONS 

In this paper we propose a probability that small glitches of pulsars could be the result of accreting SNs 
which are relics of QCD phase transition in the early Universe. This kind of glitch-trigger mechanism 
could be a possible astrophysical consequence of SNs that formed in the early Universe and survived 
to exist in our present Galaxy. The trigger mechanism is demonstrated in the framework of star-quake 
of solid strange stars. Under the typical parameters in solid strange star model, the collision between an 
SN and a pulsar could lead to the broken of the solid body (or crust) of the star, which would give rise 
to a glitch. The amount of energy releasing of such SN-triggered glitches would be smaller than that of 
the normal glitches, as the latter occur when the elastic energy accumulated reaches a critical value. In 
the star-quake model the glitch size Afi/fl is proportional to the amount of energy releasing E le , so the 
sizes of SN-triggered glitches should be smaller than that of the normal glitches. Taking into account 
that some uncertainties of both E le and the mass M for different pulsars could be mutually eliminated 
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Table 1 25 Pulsars Whose Recorded Glitch Numbers (NGlt) Are Larger Than 
5 


Pulsar name 

Spin-down age 1 
(yr) 

NGlt 2 

Distance from the 

MW center 1 (kpc) 

B0355+54 

5.64 

X 

10 5 

6 

8.87 

J0631+1036 

4.36 

X 

10 4 

15 

14.29 

B0740-28 

1.57 

X 

10 s 

8 

9.06 

B1737-30 

2.06 

X 

10 4 

33 

7.60 

B1758-23 

5.83 

X 

10 4 

12 

4.06 

J1814-1744 

8.46 

X 

10 4 

7 

2.67 

B1822-09 

2.32 

X 

10 5 

6 

7.72 

B1900+06 

1.38 

X 

10 6 

6 

5.19 

B2224+65 

1.12 

X 

10® 

5 

9.39 

J0205+6449 

5.37 

X 

10 3 

6 

10.21 

B0531+21 

1.26 

X 

10 3 

25 

9.99 

J0537-6910 

4.93 

X 

10 3 

23 

53.16 

J0729-1448 

3.52 

X 

10 4 

5 

11.30 

B0833-45 

1.13 

X 

10 4 

17 

8.04 

B1046-58 

2.03 

X 

10 4 

6 

7.66 

J1105-6107 

6.33 

X 

10 4 

5 

7.50 

B1338-62 

1.21 

X 

10 4 

23 

7.18 

J1413-6141 

1.36 

X 

10 4 

7 

6.37 

J1420-6048 

1.30 

X 

10 4 

5 

6.32 

B1757-24 

1.55 

X 

10 4 

5 

3.44 

B1800-21 

1.58 

X 

10 4 

5 

5.05 

B1823-13 

2.14 

X 

10 4 

6 

4.28 

J1841-0524 

3.02 

X 

10 4 

5 

5.50 

B1951+32 

1.07 

X 

10 5 

6 

7.46 

J2229+6114 

1.05 

X 

10 4 

6 

9.31 


( From ATNF Pulsar Catalogue jManrhpsti^ at |2QQ5l) l httn7/www atnf rsirn Ml i/ni'nnl(7ni I R:i r/rKrr:i 1 / 

“ From Jodrell Bank glitch-tables fEspinoza et alL 21) 1 lli imuf// www. ib.nian.au.uk/uuhai/aliLolicMnLini 1 


I 


to some degree, we distinguish the SN-triggered glitches from the normal glitches by their sizes with 
Afl/fl < Afl/f 2(E re = 10 35 erg), under the typical values for the mass M and radius R. 

Combined with the assumption that a huge amount of dark matter (we assume 10% in this pa¬ 
per) is in the form of SNs, and the data from the glitch data in ATNF Pulsar Catalogue and Jodrell 
Bank glitch table, we compare the collision rates of pulsars and SNs (i.e. the expect rates of glitches 
with Afl/fl < Af? /fl(£' re = 10 35 erg)) and the observational data. Among 25 pulsars that we 
choose whose recorded numbers of glitches are larger than 5, 9 pulsars have glitch events with 
An/n < An/n(£’ re = 10 35 erg) whereas 16 pulsars have no such event. The incidences of glitches 
with An/n < An/n (E le = 10 35 erg) of the 9 pulsars exhibit the tendency consistent with the colli¬ 
sion rates of pulsars and SNs, if the baryon number of each SN A ~ 10 33 . The remaining 16 pulsars 
are generally younger than the 9 pulsars, so the zero event rate might result from the fact that younger 
pulsars should release more energy during glitches. 

The mass (or baryon number A) of each SN is unknown, which depends on the properties of 
hot quark-gluon plasma and the state of dense matter. From astrophysical point of view, we want 
to give some constraints combined with observations. It seems that the constraints from Big Bang 
Nucleosynthesis (A > 10 25 ), from the supermassive black holes at high redshifts {A < 10 35 ) and 
in this paper from pulsar glitches (A ~ 10 33 ) are consistent. Certainly, SNs should have some mass 
spectrum instead of a uniform mass, so the mass corresponding to the baryon number we give above 
could be seem as the peak value of a very steep spectrum. More detailed studies about the spectrum are 
expected in the future. 
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The result we give in this paper is the first attempt to show the possibility that SNs could be a kind 
of glitch-trigger, demonstrated as the fact that there are some correlations between the glitch behaviors 
and the location of pulsars inside the Milky Way. The crude assumptions we make in this paper should 
be improved in further work, e.g. some more quantitative descriptions about the dependences of E le on 
pulsar masses and ages. In addition, a detailed description about the mechanism of small glitches as the 
result of collisions between SNs and pulsars also depends on the states of matter inside both pulsars and 
SNs. Anyhow, whether the SNs could be a kind of glitch trigger remains to be explored. 

The glitches with small sizes are now difficult for us to detect, so in some senses, to demonstrate 
the relation between SNs and glitches by the comparison between expected and observed events from 
recent data should be inadequate. However, we still try to do this because some positive results have 
been shown, which makes it possible that more data in the future would give us a clear answer. A more 
convincing demonstration about whether SNs could trigger small glitches will be done if only when 
more small glitches would be detected. This work emphasizes the importance of detecting small glitches, 
which is also the task of SKA (the Square Kilometre Array) that has exquisite timing precision and will 
dramatically increase the number of sources in surveys dWatts et all 120151) . We are then expecting to 
test our model by future advance radio facilities, e.g., SKA and Chinese FAST (Five-hundred-meter 
Aperture Spherical Telescope). 
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